COMPUTER PROGRAMS FOR KNOT TABULATION 



Charilaos Aneziris 



O 



Institut fur Hochenergienphysik, Zeuthen 
D.E.SY. (Deutsches Elektronen-SYnchrotron) 
Platanenallee 6, Zeuthen 15738 
GERMANY 



Abstract 

0\ ■ 

o\ : 

While the problem of knot classification is far from solved, it is possible to create computer 
Ch ! programs that can be used to tabulate knots up to a desired degree of complexity. Here we 
1—5 \ discuss the main ideas on which such programs can be based. We also present the actual 
00 | results obtained after running a computer program on which knots are denoted through regular 
projections. 

VO i 1. INTRODUCTION 

o 

o ■ 



In order to develop an algorithm to tabulate knots, one first needs to introduce an appropri- 
ate notation, so that knots may be presented through sequences of numbers. Such a notation 
Q\ ■ may only be valid if each acceptable presentation uniquely determines a knot up to ambient 
isotopy. One then continues by obtaining the conditions that have to be satisfied by some 
c3 | sequence of numbers, in order for the sequence to be acceptable and thus denote an actual 
qh knot. This problem is far from trivial, as shall be made clear later in the paper. 

The next step is to identify the sequences that denote ambient isotopic knots, through the 
^ . introduction of appropriate equivalence moves. If two or more number sequences are related 
through such moves, all sequences except the one appearing first are eliminated and do not 
appear at the output, and thus the corresponding knot type is tabulated exactly once. 

There is no upper bound however to the number of equivalence moves that may be needed 
to connect two sequences denoting ambient isotopic knots. Therefore if two sequences are not 
connected after some finite number of moves, this does not imply that the corresponding knots 
are inequivalent. In order to show inequivalence one needs to establish appropriate invariants. 
Each such invariant is a function on the set of notations, defined so that any sequences connected 
through equivalence moves are mapped to identical values. Therefore if one or more invariants 
map two sequences to distinct values, the corresponding knots are inequivalent. 

Ideally one may run the program long enough, so that each knot up to some desired com- 
plexity has either been shown ambient isotopic to a preceding knot, or shown inequivalent to 
all preceding knots. In the former case the knot is eliminated, while in the latter it is tabulated. 
Due however to memory and time limitations, there are certain limits up to which one may 
carry such a process. 



1 



In Section 2 we discuss two methods of denoting knots. One is based on the braid word, 
and the other on presenting knots as closed Self- Avoiding Walks (SAW's) on a cubic lattice. In 
Section 3 we present and discuss in detail the Dowker notation, which is based on studying 
knots through their regular projections. Finally in Section 4 we present the results we obtained 
after running a computer program based on the Dowker notation, and compare them to results 
obtained by others (Ref. 1). 

2. EXAMPLES OF KNOT NOTATIONS 

The Braid Word 

According to the Alexander Theorem, every knot and link may be obtained through 
the closure of appropriate braids (Ref. 2,3). Let a knot K be the closure of a braid where 
b = Y\^ =1 ap and <7j G {1, 2, n — 1}. is a generator of the braid group B n . One may thus 
denote K by the sequence ao, a±, a 2 , where ao = n and a» = 0%. Once this sequence is 
given, one may uniquely determine the braid b and its closure K. In general such sequences yield 

multicomponent links, and one thus needs a criterion to distinguish the sequences that denote 
knots from the ones that denote links. In order to do so, one replaces the braid generators oi 
with permutations pi = + 1) G Sk- The permutation P = Y\i =1 Pi belongs to a conjugacy 
class of Sk characterized by a partition p± > p 2 > ... > pk > of k. The braid closure of 
riiLiO"? is a g-component link. If q — 1, such a link is a knot. 

Two distinct sequences denote the same knot if and only if they either correspond to identical 
braids, or to braids connected through equivalence moves. The first possibility occurs by 
inserting or removing trivial factors by replacing ((XjOj^ with (CTj(7j) e for \i — j\ > 1, or 

by replacing (0-^+10^ with (cr i+l <Ji(j i+1 y . 

Two distinct braids yield ambient isotopic knots if they are connected by a series of the 
following moves (Ref. 4). 

a) b^ b'W' 1 where b, b' G B k 

(3) b <-> ba e k where b G B k , e = ±1 and ba% G B k+1 

Recently (Ref. 3) it has been shown that these two kinds of moves, known as the Markov 
moves, can be generalized and replaced by the L-moves 

b = b x b 2 <-> ab^a'afia'Y^a- 1 

where a = af...a k , a' = a^l x ...a k ^ x and b[, b' 2 are derived from bi, b 2 by replacing <x,- with <Tj +1 
Vj > %. 

Since there is no upper bound in the number of the Markov or L-moves that may be needed, 
one has to establish appropriate invariants in order to distinguish braids yielding inequivalent 
knots. Such invariants may be obtained through the skein relations 

Af()<) + BfCA = Cf{\ I) 
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from which one deduces that Af(biaib 2 ) + Bf(bia i l b 2 ) = cf{b\b 2 ). Examples of such invariants 
are the Alexander-Conway, the Jones and the homfiypt polynomials (Ref. 5). 

Such invariants however do not always distinguish inequivalent knots. In addition, knots 
expressed as braid closures may not appear in a minimum crossing number projection, as the 
knot 5-2 in Rolfsen's tables demonstrates (Ref. 6). Therefore the braid word notation becomes 
inefficient when the crossing number is large. 

Closed SAW's on a Cubic Lattice 

Let {(m,n,l)/m,n,l G Z} define a cubic lattice in R 3 . By connecting neighboring points 
on the lattice, one obtains a walk. If the initial and final points coincide, and no other point 
appears more than once, the walk is called closed Self-Avoiding, and defines a knot. A tame 
knot is defined as a knot that is ambient isotopic to a closed Self-Avoiding Walk. Therefore 
the problem of tabulating (tame) knots can be reduced to tabulating closed SAW's. 

Let ro, fl, f 2 , f* n = r~o the points on the cubic lattice defining the walk, such that \fl — 
7V. 1 1 = IVi G {1, 2, n}. The unit vectors v~l = fl — rC-i take values in {±i, ±j, ±k} which may 
be identified with the numbers {1, 2, 3, —1, —2, —3}. One may thus denote a closed SAW and 
its ambient isotopy class through a sequence of n numbers taking values in {1,2, 3,— 1, — 2,— 3} 
and indicating the vectors v~l. In order for some random sequence ai, a 2 , a n , where G 
{±1, ±2, ±3}, to denote a closed SAW, two conditions have to be satisfied. First, the sum of 
the corresponding unit vectors Y%=i ^ = in order for the walk to be closed. Second, there 
should exist no two numbers i, j, where i < j and {i,j} ^ {0,n}, such that Y?k=i+iVk — 0, in 
order for the walk to be Self- Avoiding. 

Distinct notations a iy bi yield identical walks if there is a number c G {1,2, ...,n — 1} such 
that Vi G {1, 2, n — c}, b i+c = a, and Vi G {n — c+ 1, n — 1, n}, b i+c ^ n = a i; since a switch 
from to merely relocates the end points of the walk, and is thus a reparametrization of the 
walk. 

Distinct walks yield ambient isotopic knots if they are related through the equivalence moves 
I) ^ ^ bi = a i+ i, b i+1 = ai, j ^ i, i + 1 =>- bj = aj 

II) H' ^ ^ j<i-l=>bj = dj, k = a i+ i, j >i =>• bj = a j+2 

(Possible only if ai +2 = — ai) 

While under such moves a closed SAW remains closed, it may not remain Self-Avoiding. 
One thus has to ensure that any new lattice points to be occupied after a move, were unoccupied 
before, in order for the move to be allowed (Ref. 7). 

As was the case with the braid word notation, there is no upper bound in the number of 
moves needed to connect two ambient isotopic walks, and one thus needs to establish invariants. 
In order to do so, one first performs the necessary equivalence moves so that all segments along 
one of the six directions belong to the surface of a rectangular parallelepiped that encloses the 
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walk. One then removes them and converts the closed walk to a collection of open walks, which 
are considered as the strings of a braid whose closure is the knot. It is possible from these open 
walks to obtain the braid word, and continue as before by calculating the skein invariants (Ref. 
8). 

In addition however to the problem of inequivalent knots with identical skein invariants, 
running a program based on the cubic lattice notation is extremely timing consuming. The 
simplest non-trivial knot, for example, the trefoil, needs a minimum length of 24 in order to be 
represented as a closed SAW on a cubic lattice; therefore the program will spend an enormous 
amount of time merely showing that all walks of length shorter than 24 are trivial (Ref. 9). In 
contrast, one may obtain it very fast through other notations, since its braid word is o\ G B 2 
and its crossing number is 3. 

3. THE DOWKER NOTATION 

Denoting Regular Projections 

Let a knot K and P(K) a two-dimensional projection of K. P(K) is called a regular 
projection of K if none of its multiple points has a neighborhood that looks as follows. 



I) Y II) HI) 



In other words, let f(s) be the continuous one-to-one function from S 1 to S 3 that defines 
the knot, and let t be a vector normal to the projection surface. In order for the projection to 
be normal, the following conditions must be met. 

I) Ws e S^f'is) = *&i : /'( s ) x f> 

II) V Sl ^ s 2 E S 1 : (f(s 2 ) - /00) x t = 0, (f(s 2 ) - /'(si)) x i> 

III) Sl , S 2 , S 3 e S\ Sl ^S 2 ^S 3 ^ Sl I (f( S2 ) - f( Sl j) X *1 + I (/( S3 ) - X t\ > 

In addition, in regular projections one distinguishes overcrossings from undercrossings. 

Regular projections are denoted as follows. Starting from some point and moving along the 
knot, one assigns successive natural numbers 1, 2, 2n to the crossing points according to 
the order these points are met. Each crossing is thus assigned two such numbers, ao for the 
overcrossing point and a for the undercrossing. The set of the ordered pairs (a ,a u ) denotes 
the projection. 

If one considers mirror symmetric knots equivalent and disregards connected sums, given 
an appropriate set of numbers one can uniquely determine the projection. 

Drawable — Undrawable Sets 

Let n be the crossing number. Each number in {1,2, ...,2n} appears exactly once in one 
of the n pairs. Not all such sets however denote actual projections, a simple counterexample 
being {(1,3), (2, 4)}. 
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This set, as one notices from the figure above, yields no actual projection and may thus 
be considered "undrawable" . The reason is the following. Since (1,3) belongs to the set, the 
segments 1 — 2 — 3 and 3—4 — 1 are loops. If {(1, 3), (2, 4)} were drawable, the loop 1 — 2 — 3 would 
have divided the plane into two disjoint segments according to the Jordan Curve Theorem 
(Ref. 10). Therefore the loop 3 — 4—1 would either be tangent to 1 — 2 — 3, or would have to 
intersect it at an even number of points, in order to "enter" as many times as it "exits" (vertices 
do not count). These two loops are nowhere tangent, since they share no common segment, 
and they intersect at just one point, namely at {(2,4)}. Therefore the set is undrawable. 

The general condition in order for a set to be drawable is the following. Any two loops 
defined by such a set must either share one or more segments, or intersect at an even number 
of points, vertices not counting. 

A corollary of this condition is that odd numbers must always be paired to even. If that 
were not the case, there would be some pair belonging to the set such that \j — i\ would 
be even. Therefore the loops i, i + 1, j — 1, j and j,j + 1, 2n — 1, 2n, 1, i — 1, i which 
share no common segment, would intersect at an even number of points and thus violate the 
condition. This corollary, however, while being necessary, is not sufficient, as the counterexam- 
ple {(1, 4), (3, 6), (5, 8), (7, 10), (9, 2)} shown below demonstrates. In this case, one notice that 
the loops 1 — 2 — 3 — 4 and 5 — 6 — 7 — 8 violate the drawability condition. 




In order to check whether one set is drawable, one has to consider all possible loops, which 
are 3 ra , since each crossing may or may not be a vertex of the loop. In addition, if it is a vertex, 
the loop may turn clockwise or counterclockwise. The number of pairs of loops is thus 3 2n . If 
one was to check all such pairs of loops for each possible set, the CPU time needed would be 
enormous and grow exponentially. This is not necessary, however, since the sets may be divided 
into classes of equal drawability. Sets belonging to one class are either all drawable, or all 
undrawable, and thus one has to check the pairs of loops in merely one set of each class. For 
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n = 14, for example, there are 14! (about one hundred billion) possible sets. (While checking 
drawability, one need not consider the ordering in each pair of numbers. When this is also taken 
into account, the number of possible sets becomes 2 14 = 16,384 times larger). These however 
can be divided into less than a thousand classes of equal drawability, and thus the time needed 
is significantly reduced. 

Two projections are "equally drawable" and belong to the same class, if they can be con- 
nected through the following moves. 



I) 



Add or remove a pair (i, i + and increase or decrease numbers larger or equal to i by 2. 



II) 



If a set contains pairs (i± 1, j ± 1), (i±2,j ±2) remove two consecutive pairs, decrease 
numbers between % and j by 2, numbers larger than j by 4. 



HI) 



Substitute pairs (i', k), (f, k') with pairs (i, k'), (i',f), (j, k) where \i' — — \ f — j\ 
k'-k\ - 1. 



Equivalence Moves 



According to the Reidemeister Theorem (Ref. 11), two regular knot projections belong 
to ambient isotopic knots if and only if they are connected through a series of the following 
moves. 
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I) 




These moves affect the set that denotes the projection as follows. 

First move : A pair (i,i + 1) or + is added or removed, while any number j not 
less than i is replaced by j + 2. 

Second move : Two pairs (i + 1, j + 1) or (i + 1, j), (i, j + 1) are added or removed. 

Any number k not less than mm(i,j) and less than max(ij') is replaced by k + 2, while any 
number A; not less than max(i, j) is replaced by k + 4. 

Third move : Pairs (i, j), (i', fc), (j', fc') where |i' — i\ — \f — j\ = \k' — k\ = 1, are replaced 
In (/.//).(/'./).(././,•). 

Sets connected through such moves denote ambient isotopic knots and are thus identified. 
In addition, two sets Si and S 2 are identified if 3e e { — 1,1} and c G {1,2, ...,2n} such 
that G 5*i ^> + ei mod(2n),c + e mod(2n)) e S 2 , since these two sets denote the 

same projection. The parameter c indicates a change in the origin, while e = — 1 indicates an 
orientation reversal. Finally, in order to identify mirror symmetric knots, we consider equivalent 
two sets if (i,j) G Si =>• (j, i) G S2. 



Knot Invariants 



Given a regular knot projection, one may inductively calculate any skein invariant by re- 
peatedly applying the formula 
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Af(X) + Bf(X) = Cf{\ I) 



Eventually one ends up with a sum of multiples of the invariant of the trivial knot. This 
invariant is usually set to be equal to 1. When working on the set notation however, the 
following complication arises. While switching a crossing merely changes the type of the knot, 
eliminating a crossing converts a knot into a two-component link. During the inductive process 
one obtains a number of multicomponent links, which cannot be denoted through the Dowker 
notation as defined earlier. A way to overcome this problem is to use the unordered pair {i, j} 
each time a crossing or is removed. Therefore when an unordered pair appears in 
some set, it indicates the removal of a crossing from the original knot. It is thus possible to 
keep track of the emerging links, until one ends up with unlinks consisting of trivial knots. 
Then one continues by noticing that 

/A + B\n-1 

f(n unlinked circles) = y J 

In contrast to the previous notations, one need not stop here. Inequivalent knots with identical 
skein invariants may be distinguished through "colorization" invariants. These invariants count 
the number of distinct mappings from the projection strands (segments joining undercrossings) 
to the set {1, 2, k}, such that at each crossing the three strands involved satisfy certain 
constraints (Ref. 12). The constraints are defined in such a way so that if projection becomes 
through the application of a Reidemeister move, to each acceptable mapping of Pi corresponds 
exactly one acceptable mapping of P<i- This condition is necessary in order for the number of 
mappings to be an invariant. In particular, let a strand j separate the strand i from the strand 
i + 1 as shown below, and let these strands be mapped to the "colors" q, q+i and Cj. In order 
for the mapping to be allowed, the strands must satisfy a relation Cj+i = m CiCj , where M a b is a 
k x k matrix that defines the colorization invariant. 




i+i 



In order for the number of acceptable mappings to be a true invariant, the matrix has to 
fulfill the following conditions. 

1) M„ = r V r 

2) M sr = M tr s = t 

3) M aP = 7 , M Sa = e,M yP = ri > M m = M, ■ 

Due to the enormous number of possible such matrices, it is not feasible to check them one 
after the other in order to obtain the necessary invariants. In practice we have been able to 
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obtain all matrices up to only 11 "colors", which are insufficient to distinguish knots with more 
than 7 crossings. Instead, one uses special matrices that are guaranteed to obey the above 
relations. One such category consists of matrices M a p = ta+(l — t)/3 modg where the integer 
q is prime with respect to both t and t — 1. (Corollary: q is odd) 

In order for the number of acceptable mappings to exceed the number of "colors", one 
obtains a homogeneous linear system with n — 1 equations and n — 1 unknowns, where n is the 
crossing number. This system has non-zero solutions if and only if its determinant is 0; this 
determinant is in fact the Alexander polynomial, and thus one may start by calculating this 
polynomial before actually obtaining the colorization invariants. 

A second category consists of matrices M ga9(j = gp o g a o g^ 1 , where the "colors" g p are 
elements of a conjugacy class of some group. In particular, matrices derived from permutation 
groups are extremely helpful. 

One may easily check that matrices of both these categories obey all necessary conditions 
to define invariants. 

Connected Sums 

When compiling tables of knots, one usually avoids connected sums and thus tabulates only 
prime knots. There is a number of reasons for this. First, once prime knots have been tabu- 
lated, one may obtain connected sums by combining prime knots. It is therefore unnecessary to 
consider connected sums in the computer program before the end, and one thus saves substan- 
tially in running time and computer memory. Second, connected sums of the same prime knots 
possess identical knot groups and cannot be distinguished through colorization invariants, even 
if the connected sums are distinct. Third, if one uses the Dowker notation, a set denoting the 
connected sum of k prime knots may yield up to 2 fc_1 distinct knots, and thus the notation 
becomes ambiguous. 

When one uses the braid word notation, or eventually arrives to the braid word after using 
the closed lattice notation, one may check whether a knot is prime or a connected sum as 
follows. If some braid generator <Xj e B n , where 1 < % < n — 1, appears exactly once in the 
braid word, the corresponding knot is the connected sum of the closures of the two braids that 
are formed from the products of the factors Oj for j < i and at for k > i. If both of these 
closures are non-trivial knots, the initial knot is indeed a connected sum. The problem however 
with this criterion is that if a generator appears exactly once in the braid word, the knot is not 
presented in a minimal form. Therefore one has to add at least one crossing and also increase 
the braid index by one, in order to show that a knot is a connected sum. 

This problem does not arise for the Dowker notation, which is one additional reason that 
makes such a notation preferrable. If there is some number k between 1 and 2n such that for 
any pair (i, j) belonging to the set that represents the knot, (k — 1) • (k — j) > 0, the knot is 
the connected sum of two sets, one consisting of pairs with numbers smaller than k, the other 
of pairs with numbers larger than k. If both sets denote non-trivial knots, the initial knot is a 
connected sum. One thus does not need to increase the crossing number in order to show that 
a knot is a connected sum. 
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4. COMPUTER PROGRAM RESULTS 



Using the Dowker notation we developed and ran a computer program in order to tabulate 
knots. The main ideas of the algorithm were the following. 

We first obtained all regular projections whose crossing number does not exceed a maximum 
value n. To do so, we started by considering all permutations in Sk and all fc-digit binary 
numbers, for 1 < k < n. Let S be the set denoting a projection, pi = j, hi = indicates that 
(2i — l,2j) e S, while pi = j,bi = 1 indicates that (2j,2i — 1) e S. We then eliminated the 
"undrawable" sets and the sets denoting projections that had already been recorded. 

Once the regular projections were obtained, we applied on each one of them all Reidemeister 
moves that do not increase the crossing number. If a projection was not connected to one 
checked before, the projection was stored in the computer's memory. If it was connected to 
two or more projections already stored in the memory, these projections would be identified. 
When the process was completed, the computer recorded all projections stored in the memory 
that were not identified later to any previous projections. 

When n was set 14, the number of recorded projections as a function of the crossing number 
came out as follows. (Connected sums are not included; mirror symmetric and orientation 
reversed knots are listed exactly once, even if they are not ambient isotopic) 



Number of Crossings 


Recorded projections 





1 


1 





2 





3 


1 


4 


1 


5 


2 


6 


3 


7 


7 


8 


21 


9 


49 


10 


165 


11 


552 


12 


2191 


13 


29781 



All other projections are equivalent to one or more of the recorded projections. It is also 
possible that some of the projections of the table above are ambient isotopic to each other, 
but may only be connected through Reidemeister moves involving 15 or more crossings. In 
order to obtain the distinct knot types we calculated first the Alexander polynomials of these 
projections. For projections having identical Alexander polynomials we calculated colorization 
invariants for permutation groups S p , where 1 < p < m and m is a fixed parameter. For m = 7 
we were able to show that all knots up to 11 crossings are inequivalent, while (at least) 2176 
out of the 2191 knots with crossing number 12 are inequivalent to each other and to knots with 
smaller crossing numbers (for detailed tables and explicit proofs, see Ref. 13). These results 
are consistent to results recently reported by J. Hoste, M. Thistlethwaite and J. Weeks in Ref. 
1. 
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One observes that the whole program needs merely two input numbers n and m, which 
indicate the maximum crossing number of the projections considered, and the largest permu- 
tation group used to obtain colorisation invariants. As stated earlier, we eventually set these 
parameters equal to n — 14 and m = 7. 

While running the program for smaller parameters, we also noticed the following. First, in 
order to establish all ambient isotopies for knots up to c crossings, one had to set n > y — 3 
for 5 < c < 11. To be exact, the table obtained was the following. 



c 


n 








1 





2 





3 


3 


4 


4 


5 


5 


6 


6 


7 


8 


8 


9 


9 


11 


10 


12 


11 


14 



Second, in order to prove that the tabulated knots were distinct, for c < 8 the Alexander 
polynomials were sufficient. For c = 9 or 10 one had to set m > 5, while for c = 11 one had to 
set m > 7. 
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